Chirally symmetric and confining dense matter with a diffused quark Fermi surface. 
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It is possible that at low temperatures and large density there exists a confining matter with 
restored chiral symmetry, just after the dense nuclear matter with broken chiral symmetry. Such 
a phase has sofar been studied within a confining and chirally symmetric model assuming a rigid 
quark Fermi surface. In the confining quarkyonic matter, however, near the Fermi surface quarks 
group into color-singlet baryons. Interaction between quarks leads to a diffusion of the quark Fermi 
surface. Here we study effects of such diffusion and verify that it does not destroy a possible existence 
of a confining but chirally symmetric matter at low temperatures. 

PACS numbers: 11.30.Rd, 25.75.Nq, 12. 38. Aw 



INTRODUCTION 

The QCD phase diagram is an old and very intriguing 
question. What will happen with the strongly interact- 
ing matter at large temperatures and/or densities? The 
most interesting question is the interconnection of the de- 
confinement and the chiral restoration phase transitions 
(crossovers). For many years it was considered as almost 
selfevident that the deconfmement and chiral restoration 
lines on the QCD phase diagram coincide. The reason 
for such expectation was the belief that in the confining 
mode mass of hadrons is directly related to the quark 
condensate of the vacuum. Consequently, beyond the 
chiral restoration line hadrons cannot exist and the QCD 
matter should be in a deconfmcd plasma form. At low 
densities it is indeed established on the lattice that both 
chiral restoration and deconfmement crossovers coincide 
or are rather close to each other 0, 0] • What happens 
with these lines at larger densities is unknown. 

In the large TV C world at low temperatures confinement 
persists up to arbitrary large densities [3j. This is be- 
cause both the quark - antiquark and quark - quark hole 
loops are suppressed at large TV C . Consequently, there 
is no back reaction of quarks on gluonic dynamics (no 
Debye screening of the confining gluonic field) and con- 
finement persists like in vacuum. In such system the only 
allowed excitation modes are of the color-singlet hadronic 
type. The uncorrelated single quark excitations are not 
allowed. In this case it is possible to define a quarky- 
onic matter as a dense confining matter with baryonic 
excitation modes Q. 

In the real TV C = 3 world at some large critical den- 
sity the confining gluonic field might be screened and a 
deconfining transition (crossover) would appear. Then 
a key question is how big is this critical density? Lat- 
tice simulations for the TV C = 2 QCD suggest that at 
low temperatures the deconfmement transition happens 
at densities ~ 100 times bigger than the normal nuclear 
matter density [4j . Since the TV C = 3 world is between the 
two known limiting cases (TV C = 2; TV C = oo), we expect 
that at N c = 3 the deconfmement at low temperatures 



happens at the very high densities, much larger than can 
be achieved in our laboratories or in the neuteron stars. 

Note, by definition the quarkyonic matter is a dense 
cold matter with confinement. Nothing can be apriori 
concluded about existence or nonexistence of the chiral 
restoration phase transition within such a dense matter. 
If the chiral restoration transition does exist within the 
quarkyonic matter, then it would imply that at some con- 
ditions there exists a QCD matter with confinement and 
with unbroken chiral symmetry (we imply for simplic- 
ity the chiral limit). The mass origin in such a matter 
is obviously not related to dynamical chiral symmetry 
breaking in the vacuum. 

As mentioned above, such a possibility had not even 
been considered in the past on apriori grounds. Indeed, 
the old Casher argument [j| claims that the chiral sym- 
metry breaking is required for quarks to be confined. In 
addition, it is understood that chiral symmetry break- 
ing in the vacuum is important for the mass generation 
of hadrons such as TV, p, or tt. Then, naively, hadrons 
with nonzero mass cannot exist in a world with unbroken 
chiral symmetry. 

However, the Casher argument is not general and can 
be easily bypassed |f|. Recent lattice simulations have 
convincingly demonstrated Q that in the world with- 
out the low-lying eigenmodes of the Dirac operator (i.e. 
with the artificially restored chiral symmetry) hadrons 
still exist and confinement persists. Finally, if effective 
chiral restoration in highly excited hadrons is correct [§- 



then it is possible to have hadrons with the mass 



that is not directly related to the quark condensate of 
the vacuum. 

A key question is to clarify whether existence of con- 
fining but chirally symmetric dense, cold matter is possi- 
ble. We cannot solve QCD and answer this question from 
first principles. In this situation an important step would 
be to construct a model for such a matter and see what 
mechanism could be at work. Clearly, the model must be 
manifestly confining, chirally symmetric and provide dy- 
namical breaking of chiral symmetry. Constructions that 
are based on the NJL model or the linear sigma model 
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and their extensions are not suited because of luck of 
confinement of quarks. 

The simplest possible model that satisfies all required 
criteria is the model of refs. [HI, E5| . This 3+1 dimen- 
sional model can be considered as a straightforward gen- 
eralization of the QCD in 1+1 dimension in the large N c 



limit - the 't Hooft model [1J]. It is assumed within the 



model that the only gluonic interaction between quarks is 
an instantaneous linear potential of Coulomb type. Such 
a potential is indeed observed in variational [15| as well 
as lattice Coulomb gauge simulations [l6| • An important 
aspect of this model is that it exhibits effective chiral 
restoration in hadrons with large spins J ll|, |l7|. 22 1 

This model was used in ref. [l8[ to answer the ques- 
tion about possible mechanism that could be responsible 
for the confing but chirally symmetric dense, cold mat- 
ter. Assuming a liquid phase, i.e., that the rotational 
and translational symmetries are not broken (as it is in 
the nuclear matter), the Pauli blocking of the positive 
energy levels by valence quarks prevents the gap equa- 
tion to generate a nontrivial solution with broken chiral 
symmetry above some critical quark Fermi momentum. 
At the same time confinement still persists in the system, 
so that the uncorrelated colored single quark excitations 
are not possible. 

In the latter work a rigid quark Fermi surface was as- 
sumed, like for the noninteracting fermions at T = 0. 
If such a phase exists, however, relevant degrees of free- 
dom near the Fermi surface are baryons. Quarks inter- 
act inside baryons. Consequently, the quark distribution 
function near the Fermi surface must be smooth. Here 
we study effects of such diffusion of the quark Fermi sur- 
face and solve the corresponding gap and Bethe-Salpeter 
equations. We conclude that for any reasonable diffusion 
there always exists such a critical "Fermi momentum" 
of quarks at which the chiral restoration phase transi- 
tion persists and, hence, elementary excitations above 
this critical " Fermi momentum" are the color singlet chi- 
rally symmetric hadron modes. 



OVERVIEW OF CONFINEMENT AND CHIRAL 
SYMMETRY BREAKING IN A VACUUM 

The SU(2) L x SU(2) R x U(1) A x U(l) v symmetric 
Hamiltonian 



H = 



+ - / d 3 xd 3 y J a Jx,t)K%,(x-y)4(y,t), (1) 



relies on the quark current-current J° (x, t) — 
V>(x*, i)7 M 4j- ?/>(:?, t) interaction via the instantaneous lin- 
ear interquark potential of the Coulomb type 



K£(x-if)=g fl0 g u0 5 ab V(\x-y\y, 



x a x a 



V(r) 



. (2) 

where a, b are color indices. This model was intensively 
used in the past to study chiral symmetry breakin g, ch iral 
properties of hadrons, etc, see, for example, refs. [12, LL3|, 

The Fourier transform of the linear potential and any 
loop integral are not defined in the infrared region, p ~ 0. 
Hence, an infrared regularization is required. Physical 
color singlet observables, such as hadron masses, chiral 
condensate, etc., must be independent of the infrared 
regulator /i//? in the infrared limit — > 0. 

There are several physically equivalent ways to perform 
this infrared regularization. Here we follow Ref. [20( and 
define the potential in momentum space as 



V(p) 
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(P 2 + Mi 2 r) 5 



(3) 



This potential in the configuration space contains the 
required or term, the infrared divergent term —cr/^in. a $ 
well as terms that vanish in the infrared limit. 
The self-energy operator 



(4) 



consists of the Lorentz-scalar chiral symmetry breaking 
part A p and chirally symmetric part (jp) (B p — p) . The 
unknown functions A p and B p are to be determined from 
the gap equation for the chiral angle ip p 



A p cos ip p — Bp sin tp p = 0, 



where 



A„ 



1 f d 3 k 



2 J (2tt) : 



V(p- k) sin^, 



(5) 



(6) 



1 f d 3 k 

B p =p+ 2 J (2tt)3 (P k ^ V (P~ k)costp k . (7) 

These integrals contain both the infrared divergent and 
the infrared finite parts 
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The same is true for the single quark energy: 



(8) 
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(10) 



Consequently, the single quark Green function is diver- 
gent and the single quark energy is infinite in the infrared 
limit. Actually, energy of any color nonsinglet state is 
infinite. At the same time the infrared divergence ex- 
actly cancels out in any color singlet quantity and these 
quantities are finite and well defined 18J . This is a man- 
ifestation of confinement within this model. 

Similarly, the infrared divergence in A p and B p can- 
cels in the gap equation and this equation can be solved 
directly in the infrared limit. The gap equation can be 
solved numerically and a nontrivial solution with the chi- 
ral (Bogoliubov) angle ip p ^ signals dynamical breaking 
of chiral symmetry. Hence, the single quark Green func- 
tion is not chirally symmetric, A p ^ 0; there appears the 
quark condensate and dynamical mass of quarks 
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CHIRAL SYMMETRY BREAKING AND 
CONFINEMENT IN A DENSE MATTER AT T=0 
WITH A RIGID QUARK FERMI SURFACE 

It is practically impossible to solve exactly the model 
in a dense matter. Indeed, that would imply to solve it 
first for a single baryon; then to obtain a baryon-baryon 
interaction; given this interaction to construct a nuclear 
matter and then slowly to increase its density. Obviously, 
it is a formidable problem. In order to proceed and get 
some insight one needs justifiable simplifications. 

In the large N c limit the nucleon is infinitely heavy, 
translational invariance is broken and a many-nucleon 
system is certainly in a crystal phase. Whether a (dense) 
nuclear matter will be a liquid or a crystal at N c = 3 is 
a subject to dynamical calculations. Such microscopical 
calculations cannot be persued for any "realistic" model 
in 3+1 dimensions with confinement and (broken) chiral 
symmetry. However, in the real world iV c = 3 we do 
know that the nuclear matter is in a liquid phase; both 
translational and rotational invariances are intact. We 
then assume a liquid phase with manifest translational 
and rotational invariances in a dense quarkyonic matter. 

We also assume that confinement persists up to large 
densities at N c — 3. At T = deconfinement could 
happen through the Debye screening of the confining 
gluon propagator: A gluon creates the quark - quark 
hole pair that again annihilates into a gluon. If this vac- 
uum polarization diagram is finite, then at some density 
there should happen a complete screening of the confin- 
ing gluon field. However, in the confining mode energy 



of the colored quark - quark hole pair is infinite |2l| . The 
allowed excitations in the confining mode are the color 
singlet excitations like baryon - baryon hole pairs, etc. 
These excitation modes cannot screen the confining col- 
ored gluon propagator. In this sense the T = physics 
is rather different from the deconfinement at zero den- 
sity and large temperature. In the latter case a screening 
proceeds via the incoherent thermal gluon loops. 

One could expect that the deconfinement should hap- 
pen in a dense medium at T = due to perlocation of 
baryons. Such a reasoning is too naive, however, be- 
cause the perlocation does not yet imply screening of the 
confining gluon field. For example, deconfinement never 
happens at T = in QCD at large N c or in the 't Hooft 
model. In both cases baryons "sit" on top of each other 
in a very dense medium, but it is still a system with 
confinement. 

We want to address the chiral symmetry breaking 
properties of a dense matter with confinement. In the 
vacuum dynamical chiral symmetry breaking happens 
because there is an attractive interaction between the left 
quarks and the right antiquarks and vice versa. This at- 
tractive interaction shifts the energy of the vacuum with 
broken chiral symmetry below the energy of the pertur- 
bative Dirac vacuum. 

Consequently, in a dense matter at T = the most 
important physics that leads to the restoration of chiral 
symmetry is the Pauli blocking (by the valence quarks) of 
the positive energy levels required for the very existence 
of the quark condensate. 
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FIG. 1: Valence quark distribution and the corresponding 
integration weight. 

In our previous work [l8[ the effect of the Pauli block- 
ing was studied assuming a rigid valence quark Fermi 
sphere, like for the ideal Fermi gas at T = 0, see Fig. 
la. In this case all quark positive energy levels below 
the Fermi momentum pf are occupied by the valence 
quarks and one has to replace the vacuum density matrix 
v(p)v^(p) by the density matrix in the medium: 



Pip) = ®(Pf ~ p)u{p)v) (p) +v(p)v\p). 



(12) 



Hence, within the mean field approximation one has to 
remove from the integration in the gap equation all quark 
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FIG. 2: Quark condensate in units of a 3 ^ 2 as a function of 
the Fermi momentum, which is units of y/a. 

momenta below pf since they are Pauli blocked, see Fig. 
lb. The modified gap equation is then the same as in ([5]) 
- ([7]). but the integration starts not from k = 0, but from 
k=pf. 

At the critical Fermi momentum p c ^ ~ 0.109-y/cr the 
chiral phase transition is observed, see Fig. 2, because 
the nontrivial solution with broken chiral symmetry dis- 
appears. Above this phase transition the chiral angle ip p , 
the quark condensate (qq) , the dynamical mass of quarks 
M(j>) as well as the chiral symmetry breaking part A p of 
the quark Green function identically vanish. At the same 
time the chirally symmetric part B p of the quark Green 
function does not vanish and is still in fact infrared diver- 
gent. The single quark energy is infinite and confinement 
persists even in the chirally symmetric phase. The color 
singlet hadronic excitations have finite and well defined 
energy. 

EFFECTS OF A DIFFUSION OF THE QUARK 
FERMI SURFACE 




a) b) 

FIG. 3: Diffused step valence quark distribution and the cor- 
responding integration weight. 

In reality valence quarks near the Fermi surface inter- 
act and cluster into the color singlet baryons. This in- 



teraction in general would lead to a diffusion of the rigid 
Fermi surface for quarks. Some levels above the "Fermi 
momentum" must be occupied with some probability as 
well as some levels below the "Fermi momentum" with 
some probability must be empty. 

In principle, the quark distribution function near the 
diffused Fermi surface could be obtained selfconsistently 
from the full solution of the problem. It is a formidable 
task and such a program cannot be persued. However, 
it is clear that the realistic distribution function will be 
smooth, of the form on Fig. 3a. Then, for our present 
goal it will suffice to parameterize such a smooth distri- 
bution in a simple form and study effect of the diffusion 
on the solution of the gap equation. 

We parameterize a smooth valence quark distribution 
function by 

p p (p) = Q(- P +p f -A) + e(p -p f + A) eiP _J /A - 1 ■ 

(13) 

Given the valence distribution function we multiply the 
integrands in eqs. © - © by the weight function 1 — 
p v (k) and solve the gap equation for different pf and 
diffusion width A. 

If the diffusion width is much smaller than the crit- 
ical Fermi momentum for a rigid Fermi surface, A <C 
p c f(A = 0), (what should be considered as a realistic sit- 
uation) , then the evolution of the chiral condensate with 
Pf is similar to the case of the rigid quark Fermi sur- 
face. This situation is represented by the curve A = 0.02 
on Fig. 4. The phase transition happens at the "Fermi 
momenta" that are rather close to the critical Fermi mo- 
mentum, p C f{A = 0) = 0.109 from Fig. 2 (see also the 
curve A = on Fig. 4). This can be easily understood. 
At all momenta p <C p/ the Pauli blocking on Fig. 3 is 
the same as for the rigid quark Fermi surface. At mo- 
menta just below the pf the effect of the Pauli blocking 
is weaker than for the rigid Fermi surface. However, this 
is compensated by additional Pauli blocking of the levels 
that are just above the Fermi momentum for the rigid 
quark distribution. Consequently, with a small diffusion 
widths the "critical Fermi momentum", p^ r (A), at which 
the phase transition happens, is shifted to slightly lower 
values oi pf. 

However, if a diffusion width becomes larger and even- 
tually comparable with the " Fermi momentum" , then the 
" critical Fermi momentum" , at which the phase transi- 
tion happens, increases. 

For each fixed diffusion width A there always exist 
such "Fermi momenta" where the Wigner-Weyl mode of 
chiral symmetry is realized. This can be seen from Fig. 5, 
where a line of " critical Fermi momenta" is depicted. The 
area above this critical line corresponds to the chirally 
symmetric phase, while all points below the critical line 
represent a matter with broken chiral symmetry. 

The chiral angle <p p and dynamical mass of quarks 
M(p) in the Nambu - Goldstone mode of chiral symme- 
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FIG. 4: Quark condensate (in units of <r 3 / 2 ) as a function of 
the Fermi momentum (in units of y/cf) and the diffusion width 
A (in units of \fa) 



FIG. 6: Chiral angle tp p as a function of the momentum p 
for a "Fermi momentum" pj — 0.2 at different fixed values of 
smoothing A. p and A are in units of yfa. 




FIG. 5: Critical line that separates the quarkyonic matter 
with broken and restored chiral symmetry. 



try at some "Fermi momentum" and different diffusion 
widths A are shown on Figs. 6-7. 



CHIRAL RESTORATION IN MESON SPECTRA 

Another explicit illustration of chiral symmetry of a 
dense matter above the chiral restoration phase tran- 
sition are properties of hadronic excitations. In the 
Nambu-Goldstone mode of chiral symmetry there must 
be a massless excitation mode that is associated with the 
massless pion. At the same time energies of all other 
mesons must be finite. In particular, there must be a 
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FIG. 7: Dynamical mass M(p) as a function of the momentum 
p for a "Fermi momentum" p/ = 0.2 at different fixed values 
of smoothing A. M(p), p and A are in units of yfa. 



finite splitting of the excitations with quantum numbers 
I, J pc = 1,0"+ and I, J PC = 0,0++, that will be ref- 
fered as the pion and the er-meson, respectively, according 
to the standard nomenclature. In contrast, these exci- 
tations must be exactly degenerate in the Wigner-Weyl 
mode of chiral symmetry and form the (1/2, l/2)„, repre- 
sentation of the SU(2) L x SU(2) R chiral group [HI ]. 

To obtain the quark- antiquark bound states we solve 
the homogeneous Bethe-Salpeter equation in the rest 
frame 



X(m,p) 



d 4 q 
(2tt)4 



V(\p-q\) 7o S(q + m/2,p-q) 
P°(3)).(14) 
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Here S is the dressed single quark propagator, that is 
the solution of the gap equation, m is the meson mass 
and p is the relative momentum. The Bethe-Salpeter 
equation is solved by means of expansion of the vertex 
function x{ m iP) m to a set of all possible independent 
amplitudes consistent with /, J PC and it transforms into 
a system of coupled equations. The infrared divergence 
cancels exactly in these equations and they can be solved 
numerically [17| . The Pauli blocking by valence quarks 
is taken into account via the weight function 1 — p v (q) . 

In the Wigner-Weyl mode, i.e., when dynamical quark 
mass and chiral angle vanish, M(j>) — 0; ip p — 0, the 
Bethe-Salpeter equations for the 1, h and 0, ++ bound 
states become identical [T3| and consequently energies of 
these states coincide. 

On Fig. 8 we show masses of both pseudoscalar and 
scalar modes for different "Fermi momenta" pf and dif- 
fusion widths A. For each A there is a critical pj r (A) at 
which the chiral restoration phase transition takes place. 
Below this pj r (A) there is a massless pion and a massive 
cr-meson. Above the critical pj both the pion and the 
tr-meson are massive and exactly degenerate. 



CONCLUSIONS 

In the confining mode the valence quarks interact 
and near the Fermi surface cluster into the color singlet 
baryons. This implies that there cannot be a rigid quark 
Fermi surface. The valence quark distribution function 
near the Fermi surface must be smooth. The valence 
quark levels above the "Fermi momentum" are occu- 
pied with some probability as well as the levels below 
the "Fermi momentum" must be, with some probability, 
empty. We assume unbroken translational and rotational 
invariances, i. e., a liquid phase. We parameterize such 
a diffused "Fermi surface" by a simplest possible func- 
tion and solve the corresponding gap and Bethe-Salpeter 
equations. By this we verify whether a chiral phase tran- 
sition, previously observed for a rigid quark Fermi sur- 
face, survives or not. It turns out that for any reason- 
able diffusion width there always exists such a "Fermi 
momentum" that the chiral restoration phase transition 
does take place. This reconfirms our previous conclu- 
sions about possible existence of the confining but chi- 
rally symmetric phase. Below the phase transition the el- 
ementary excitation modes of a matter are hadrons with 
broken chiral symmetry, while above the phase transition 
such excitations are chirally symmetric hadrons. 
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FIG. 8: Masses of the pseudoscalar (solid) and scalar (dashed) 
mesons in units of \fa as functions of the "Fermi momentum" 
Pf and of the diffusion width A (in units of 
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